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Abstract. The evidence for dual superconductivity as a mechanism for color confine- 
ment is reviewed. New developments are presented for full QCD, i.e. in the presence of 
dynamical quarks. 



1. Introduction 

Confinement is by definition the absence of free colored particles in nature. 
In spite of the clear evidence coming from high energy experiments that 
quarks and gluons are the fundamental constituents of hadrons, none of 
them has ever been detected. 

Only upper limits exist to production cross sections .The cross section 
aq for the inclusive process 

p + p^q{q)+X (1) 

has an upper limit 

ag < 10-4° cm2 (2) 
At the same energies the total cross section or has the value 

ar ^ lO^^^cm^ (3) 

In perturbation theory the ratio Cq/aT is expected to be a sizable fraction 
of unity. From (2) and (3) 

aq/ar < IQ-^^ (4) 

Relic quarks in nature have been hunted since they were first proposed as 
fundamental bricks of matter [01. 
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Fourty years of Millikan-like experiments looking for fractionally charged 
particles have produced as upper limit [1] 

n^/rip < 10^27 

for the ratio of the abundance of quarks Ug and that of nucleons Up in 
nature. The limit (5) results from the analysis of ~ Ig of matter and no 
quarks found. In the absence of confinement the Standard Cosmological 
Model predicts 

ng/up ~ 10^12 (g) 

Again 

{nq/np)obs ^ io~i5 (7) 



{jlq / TLp^ expected 

The only natural explanation of the limits (3) and (7) is that cjg, Ug are ex- 
actly zero, or that confinement is an absolute property based on symmetry 
[4]. Similar situations are e.g. 

i) the photon mass, which is experimentally bounded by the inverse 
radius of the solar system[p.The corresponding symmetry is gauge 
invar iance. 

ii) ordinary superconductivity. The upper limit of the resistivity of super- 
conductors is many orders of magnitude smaller than that of any other 
material. The symmetry pattern behind that is the Higgs phenomenon^] 



2. Virtual tests (lattice) 



QCD at finite temperature can be simulated on the lattice. It is a well 
known theorem that the partition function of a field system is equal to the 
euclidean Feynman integral, with immagianary time ranging from to l/T, 
and periodic boundary conditions for bosons, antiperiodic for fermions: 



Z = i \P'^\ exp 



l/T j- 

dt I cfixC{x,t) 



If the lattice size in the time direction is Nt then 

1 



Nra 



(8) 



(9) 



a being the lattice spacing in physical units. Renormalization group gives, at 
large enough f3 = 2N/g'^, a = exp(— /3/6o), where 6o > by asymptotic 
freedom, i.e. 

T=^eMP/bo) (10) 
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Low T corresponds to large g'^ (strong coupling or disorder in the language 
of statistical mechanics ); high T to weak coupling or order. This is the 
opposite to what happens in ordinary spin systems, for which T plays the 
role of . 

In pure gauge theories an order parameter for confinement is {L), the 
Polyakov line, which is the trace of the parallel transport along the immag- 
inary time axis from to 1/T, closed by periodic boundary conditions. The 
corresponding symmetry is Z^- 

On the lattice the correlator 

G(f) = (L(f)Lt(0)) (11) 
is measured |6|. Cluster property requires 

G{f) ^^^^Aexpi-aaNxr) + \{L)\'^ (12) 

The potential energy of a static qq pair at distance r is given by 

V{r) = --^lnG{r) (13) 

A temperature Tc is found such that for T <Tc {L) = or 

V(r) ~ ar (14) 

which means confinement. For T > Tc {L) ^ and 

V{r) ~ const. (15) 



which means deconfinement. 

Finite size scaling analysis of the correlator around the critical point 
provides a determination of the critical index z^. For SU (2) pure gauge 
theory the transition is second order |^, consistent with the class of univer- 
sality of the 3d ising model {v = .62) as expectedQ, and T^j ^ ^ .7. For 
S'C/(3) pure gauge theory the transition is weak first order ^, {y = -33) 
and Tel \fo ~ .65, which, by the usual assumption = 425 MeV gives 
Tc ~ 270 Mev. 

In the presence of dynamical quarks is explicitely broken , and {L) 
cannot be an order parameter. For two equal-mass dynamical quarks the 
situation is depicted in fig.l. 

The transition temperature is determined, at given quark mass, by look- 
ing at the maximum of a number of suceptibilities , e.g. / {'ipip{x) ipip{0))d^x, 
J {L{x)L(0))d^x. All of them show a maximum at the same Tc|jl^. For high 
enough niq, {niq > 3 GeV) the maximum of the Polyakov line susceptibility 
goes large with increasing spatial volume as in the quenched case: a finite 
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Figure 1. Phase dyagram of QCD:shaded region is confined. 



size scaling analysis shows that the transition is first order. There are in- 
dications that the transition is second order in the chiral limit niq = 0,as 



suggested by symmetry arguments[ll|. At intermediate values of niq none 



of the susceptibilities which have been considered increases with increasing 
spatial volume, and a possible conclusion is that there is no phase transition 
but only a crossover. The overall situation is rather confusing. It is not clear 
a priori what is the relation between chiral symmetry and confinement. It 
is not fully clear either what susceptibilities are entitled to determine the 
order of the transition by their behavior at large volumes. In principle the 
relevant quantities should be those appearing in the expression of the free 
energy. The free energy (effective lagrangean ) depends on the dominant 
exitations and on their symmetry. What are the dominant exitations is 
exactly the problem under investigation. 

3. Duality. 

Confined phase is disordered. How can the symmetry of a disordered phase 



be defined? The key concept is duality [ 12 1. It applies to d-dimensional 
systems admitting non trivial topological excitations in (d-1) dimensions. 
These systems admit two complementary descriptions. 

1) A direct description in terms of the fields </>,with order parameters (0), 
in which the topological configurations /i are non local. This description 
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is convenient in the weak coupling regime {g <ti 1), i.e. in the ordered 
phase. 

2) A dual description in which the topological excitations /i become local 
fields, and the original fields (p topological configurations. The dual 
coupling gjj is related to g as go ^ 1/5- This description is convenient 
in the disordered phase (strong coupling regime). Its symmetry is de- 
scribed by (/x) (disorder parameter). Duality maps the strong coupling 
regime of the direct description into the weak coupling regime of the 
dual description. 

The prototype system for duality is the ising mo del [ p^ w here dual excita- 
tions are kinks. Other examples are N=2 SUSY QCD|14||, where the dual 
excitations are monopoles; M string theories fisll ; 3-d XY model, where dual 
excitations are abelian vortic es [[iq ]; 3-d Heisenmerg magnet, with 2-d Weiss 
domains as dual excitations IT^I; compact U(l) gauge theory , where dual 



excitations are monopoles [ 15 , 19]. 



In QCD the dual topological excitations have to be identified : as we will 
see, however, information exists on their symmetry. Two original proposals 
exist in the literature, which have been widely studied: 



a) Monopoles ]]2^, 21]. The idea is that vacuum acts as a dual supercon- 
ductor, which confines electric charges by Meissner effect, in the same 
way as magnetic charges are confined in an ordinary superconductor. 
Developments of this approch will be the subject of the next sections. 

b) VorticesQ. The symmetry involved is Z^. 

In 2-|-l dimensions a conserved charge exists, the number of vortices minus 
the number of antivortices, and vortives are described by a local field. In 
3-1-1 dimensions a dual Wilson loop can be defined ('tHooft loop ) B{C), in 
connection with any closed path C . The algebra which is obeyed by B{C) 
and by the ordinary Wilson Loop W{C') is 

B{C)W{C') = W{C')B{C) exp (incc ^"l (16) 



where ncc is the linking number of the two loops. From eq(16) it follows 
that, if {W{C')) obeys the area law {B{C)) obeys the perimeter law, and 
if {B{C)) obeys the area law then {W{C')) obeys the perimeter law. If we 
denote by {L) the ordinary Wilson loop which wraps the lattice through 
periodic b.c. in time (Polyakov loop) ,and by (L) the analogous dual loop 
('tHooft's line), then in the confined phase {L) = 0, (L) ^ 0, whilst in 
the deconfined phase (L) = 0, (L) / 0. (L) is a disorder parameter for 
confinement. These relations have been tested on the lattice[^, The 
corresponding symmetries and Z^ are explicitely broken in the presence 
of fermions. 
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4. Monopoles. 

Monopoles in non abelian gauge theories are always abelian (Dirac) monopoles. 
This statement can be immediately checked by looking at the field produced 
by a static configuration of colored matter at large distances, by use of the 
familiar multipole expansion [24]. Monopoles are identified by a constant 
diagonal matrix in the algebra, with integer or half-integer values : they 
carry N-1 abelian magnetic charges. The same physics emerges from the 
procedure known as abelian projection [21]. We shall illustrate it for SU{2): 
the general case|2^] is not substantially different. Let (p{x) be any operator 
in the adjoint representation, and (^(x) = ip{x)/\(p{x)\ its direction in color 
space. Define [p6[ 

- 1 

i^M^ = ^G^^ - -^{Df^ip A D^if) (17) 

g 

with Gfj_i, = d^Aj^ — dyA^ + gA^f\Ay the field strength and = d^ + gA^f\ 
the covariant derivative. Both terms in eq.(17) are color singlets and gauge 
invariant: the combination is chosen to cancel bilinear terms A^A^. Indeed 
one has identically: 

F^,u = vid^Ajy - dyA^) - ^ifid^if A di^if) (18) 

In a gauge in which is constant, e.g. = (0,0, 1), F^i, is abelian : 

Ffxu — d^Aj^ — duA^ 

A magnetic current , j^, can be defined in terms of the dual tensor F*^^ = 
ie F 

2^^i'p<y^ per? 

3^l = dyF^^ 

is identically zero ( Bianchi identities ) in a non compact formulation of 
the theory. In a compact formulation, like Lattice, can be non zero. In 
any case it is identically conserved 

dt^f = (19) 

Magnetic charges are Dirac monopoles, obeying Dirac quantization condi- 
tion Q = n/2g. The corresponding magnetic U{\) symmetry can either be 
realized a la Wigner, and then the Hilbert space consists of superselected 
sectors with definite magnetic charge, or Higgs-broken, and then the system 
behaves as a dual superconductor. If the ideas of ref 's | pO| , p| ] are correct the 
expectation is that QCD vacuum behaves as a dual superconductor (Higgs- 
broken phase) for T < Tc, and as a magnetic superselected system for 
T > Tc- A disorder parameter should discriminate between superconductor 
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and normal. Such a parameter has been constructed! 27, 28, ^] as the v.e.v. 
{jj) of an operator ^ carrying magnetic charge. In fact is a Dirac-hke 
operator [30 1, charged and gauge invariant [p^, |33 ]. The construction of n is 
at the level of a theorem for compact U (1)|18[ |33|]. In non abelian gauge 
theories it is undefined by terms O{o?'), a being the lattice spacing, like the 
abelian projection itself 1 32, p^ ]. 

5. Results 



The basic structure of /i is a translation of the field configuration in the 
Schrodinger picture by a classical monopole configuration. In the same way 
as 



e^P"-\q) = \q + a) 



(20) 



defining 



^(x, t) =exp[i d yn{y, t)ip{x - y) 



with n(y, t) the conjugate momentum to the field f{y,t), and ip{y,t) the 
classical field configuration to be added 



IJ.{x,t)\ip{y,t)) = \ip{y,t) + ip{x-y)) 



(21) 



In fact the basic structure has to be adapted to a compact formulation, in 
which the field cannot be translated at will|jl9[, and to the abelian projected 
situation, in which only the abelian part of the field has to be translated. 



All this has been done[27|. The resulting disorder parameter (/i) can be 
finally written as the ratio of two partition functions 



m 



(22) 



with Z(/3) = J[Vip] exp(-/35), and 

ZiP) = J [V^] exp [-(3{S + AS)] (23) 

S + AS" is obtained from S by a modification of the space time plaquettes 
at time t, IlQi{n,t); for details see ref's |19, Instead of (/i) itself it is 
more convenient to study 

d 



ln{fi) = {S)s-{S + AS)sAS 
dp 



(24) 



On the one hand p is numerically easier, since {fi) fluctuates wildly as any 
partition function; on the other hand we shall see that p contains all the 
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relevant information as (/i) , and also in a more convenient way. From p, (fi) 
is obtained as 

r/3 



(p) = exp 



p{x)dx 



(25) 



since Z{f3 = 0) = = 0) = 1. 

The typical shapes of {p) and p as functions of /? are plotted in fig.2. 
The position of the negative peak coincides with the deconfining phase 
transition. 




Figure 2. Typical shape of and rho. 



Fig. 3 shows p for different spatial sizes Ng of the lattice, at fixed Nt = 4, 
for SU{2) pure gauge theory. 

The position of the peak coincides with the maximum of the suscep- 
tibility of the Polyakov line, as determined in ref|^, i.e. with the phase 
transition. In the range of temperatures T < Tc p stays practically constant 
by increasing the volume, as shown in fig. 4 ,and this means that (p) has a 
non zero limit in that region. 

In the range T > p diverges to — cxd as A'^s goes large, as 



-kNs + k' 



k>0 



(26) 



as shown in fig. 5 . 

If we had measured {p) itself instead of p we would have found it 
consistent with zero within large errors at large Ns- Measuring instead 
p and checking numerically the behavior eq.(31) amounts to state that (p) 
is exactly zero in the thermodynamical limit. In the critical region, around 
Tc one expects 



(27) 



with T = 1 — ^ oc {Pc — (3), S a critical index, a the lattice spacing and ^ the 
correlation length. The transition is known to be second order for SU{2), 



9 



200 



100 



-100 



^ -200 



-300 



-400 



-500 



-600 









A D a ■ A 












m ^ 


• 


• 


V. 






? 




i N,=12 






• N=16 






T N =20 






■ N=24 


1 




*N=32 









Figure 3. 5(7(2): p at different spatial sizes. 



weak first order for 5f7(3) gauge theory, and therefore ^ ~ r goes large 

at small r's. Neglecting a/^ ~ and trading the variable Ng/^, with tN}^'^ , 
the scaling law follows for p from eq.(32) 



(28) 



For different sizes of the lattice p/n}^'^ 
expected to be independent of A'^g. A best fit to the data||2 
a determination of /3c, i^, 5. The quality of the scaling is shown in fig. 6 for 
SU{2). The result is 



when plotted versus tNs' is 

1 allows 



28 



SU{2) 
SU{3) 



.62(1) 
.33(1) 



5 
6 



.20(3) 
.50(3) 



All that is for quenched theory. A few different abelian projections 
have been tested, and the behavior of p ,as well as the value of and 
5 are independent of the abelian projection. An additional test has been 
made with no abelian projection, i.e. by assuming as diagonal operators 
the nominal A3, Ag used in the simulation!^]. This amounts to perform an 
average over all abelian projections, and the result does not change. 

We can thus state that the confining phase of a pure gauge theory 
is a dual superconductor in all the abelian projections, and undergoes a 
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Figure 5. p at T > Tc as a function of spatial size. 



transition to normal at Tc- What we learn from that is that, whatever 
the dual excitations are, they carry magnetic charge in all the abelian 
projections. The definition of the operator // can be easily extended to full 
QCD, with dynamical quarks. A natural question is whether also in this 
case the confining phase is characterized by dual superconductivity. This 
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Figure 6. Finite size scaling of p 

would indeed be the expectation in the spirit of the ^ oc hmit. If the 
physics of gauge theories is determined by the hmit Nc ^ co aX g^N^ fixed, 
corrections 1/Nc being smaU, then the mechanism of confinement has to be 
Nc independent, and also insensitive to the presence of dynamical quarks, 
their contribution being non leading in the l/Nc expansion. The parameter 
p = A In(^) should go to — oo for T > Tc in the thermodynamical limit 
so that (/i) = 0: this is indeed the case |Q. For T < Tc p converges to 
a finite limit, i.e. (/x) ^ and there is superconductivity ||3^. The shaded 
area of fig.l does indeed correspond to a superconductor, the upper area 
to a superselected magnetic system: the negative peak is at the transition 
as defined in sect. 2 [fig. 7] 

Around Tc a finite size scaling analysis can give information on the order 
of the transition : the difference with respect to the quenched case is that 
now an extra dimensionful quantity, the quark mass, enters and there is a 
two variable scaling. As usual we expect 

if')r^^^'K^^^^'^Ns) (29) 

In the critical regime a/^ ~ , Nt/Ns — 0, ^ ~ and we can again 
trade ^/Ng with tN}^'^ , and the scaling law becomes 

{p)^^^r'm^,rn,N^) (30) 



The index a is known: one can chose for different values of Ng suitable 
quark masses, so as to keep mgN" fixed. The scaling law is then the same 
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Figure 7. 2 flavour QCD. Open circles represent the chiral condensate, with the scale 
axis on the left. Pull circles represent p. The peak of p coincides with the drop of (V'V')- 



as cq (28) , whence the index i' can be extracted and with it information on 
the order of the transition. This is a heavy numerical program which is on 
the way on our APE machines. The result will possibly be relevant to the 
problems discussed in sect 2. 



6. Conclusions and outlook. 



1) Confinement is characterized by dual superconductivity of the vacuum 
in all the abelian projections, both in quenched and in full QCD, in 
line with the idea of the limit Nc — ^ oo. A universal disorder parameter 

has been defined for it. 

2) We do not know the dual excitations of QCD , nor their effective 
lagrangean. However we know that they carry magnetic charge in all 
the abelian projections. 

3) An analysis of the critical region in full QCD is on the way , which will 
possibly clarify the nature of the phase transition depicted in fig.l. 

4) The determination of the parameters of the dual supercondutor , Higgs 
mass, penetration depht,. . . is fundamental and relevant to understand 
the deconfinement signsls which could come from heavy ion collisions: 
this is also part of our research program. 

The contributions of L. Del Debbio, M. D'Elia , B. Lucini, G. Paffuti to 
the research program are aknowledged. 
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